Mon. Not. R. Astron. Soc. OOP. [TlfTTl C2007') 



Printed 2 February 2008 



(MN I*TeX style file v2.2) 



Optical emission line properties of Narrow Line Seyfert Is 
and comparison AGN 

J. R. Mullaneyi* and M. J. Ward^ 

^ Institute for Computational Cosmology, Durham University, South Road, Durham, DHl 3LE, U.K. 
Date Accepted 

ABSTRACT 

Based on a new spectroscopic sample observed using the WHT, we examine the kine- 
matic properties of the various emission hne regions in narrow hue Seyfert 1 galaxies 
(NLSls) by modelling their profiles using multiple component fits. We interpret these 
results by comparison with velocity components observed for different lines species 
covered in the same spectrum, and equivalent components measured in the spectra 
of some broad line Seyfert Is and a representative Seyfert 2 galaxy. We find that the 
fits to the Ha and H/3 line profiles in NLSls require an additional broad (>3000km/s) 
component that might correspond to a suppressed broad line region with similar kine- 
matics to those of typical broad line Seyfert Is. From the profiles of the forbidden 
high ionisation lines (FHILs) in NLSls, we find evidence that they appear to trace an 
'intermediate' velocity region with kinematics between the standard broad and narrow 
line regions. Weaker evidence of this region is also present in the profiles of the per- 
mitted Balmer lines. Finally, we note that despite having similar ionisation potentials, 
the relative intensities of the highly ionised lines of [Fe XJA6374 and [FeXI]A7892 show 
considerable dispersion from one galaxy to another. The interpretation of this requires 
further modelling, but suggests the possibility of using the ratio as a diagnostic to con- 
strain the physical conditions of the FHIL emitting region and possibly the shape of 
the spectral energy distribution in the vicinity of 200eV. This spectral region is very 
difficult to observe directly due to photoelectric absorption both in our Galaxy and 
intrinsic to the source. 
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1 INTRODUCTION 

Narrow line Seyfert 1 (NLSl) galaxies are a sub-set of the 
AGN popu lation which share a number of characteristic 
properties l|Pogge|[2000l ). Superficially, their optical spectra 
resemble those of type 2 Seyferts in that their permitted 
emission lines are quite narrow and only marginally wider 
than their forbidden lines - a property that has unfortu- 
nately led to nu merous misidentifications as Sy2s in the liter- 
ature. However, lOsterbrock fc Poggel (|l985l ) noted that their 
strong permitted Fe II lines and high H/3/[0 IIIJA5007 ratios 
were similar to the traditional type 1 Seyferts. Although the 
velocity widths of the permitted lines in AGN spectra are 
observed t o cove r a continuous range (~10^ to ~10* km/s), 
iGoodrichI (|l989l ) placed a somewhat arbitrary limit of to 
the FWHM of H/3 of 2000km/s, for an object to be classed 
as a NLSl. Since then most studies have continued to use 
this value as the boundary between NLSls and their broad 
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line counterparts, here referred to as broad line Seyferts, or 
BLSls, in order to distinguish them from NLSls. 

Studies of the continuum emission from NLSl revealed 
interesting features across the electromagnetic spectrum, for 
example, significant soft X-ray excesses, rapid X-ray vari- 
ability and strong emission from high ionisation forbidden 
lines (FHIL s) re.g. lStephenslll989l : |Puchnarewicz eraLlll992l : 
|Pogge|[2000l 'l. Their X-ray properties encouraged searches us- 
ing X-ray surveys with follow-up confirmation from optical 
spectroscopy. This showed that around 50% of all ROSAT 
soft X-ray selected AGN were NLSls (|Grupe et al.l Il994l : 
iPoggdlioOOl ). It should be noted, howeve r, that this method 
introduces serious selection efltects, since iBoller et all (Il996l ') 
found that not all NLSls display strong soft X-ray excesses. 

The features described above, in particular the rapid 
X-ray variability and lack of significantly broadened per- 
mitted emission lines has led to speculation that NLSls 
may harbour a lower mass central black hole than is typ- 
ical in BLSls (1O'^ ~^M0 and lO^'^M©, respectively; e.g. 
iMathur et al.l 1200 ll ). This, coupled with their bolometric 
luminosities, implies a high accretion rate, possibly super- 
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Eddington suggesting that the key parameter in determining 
the distinctive physical properties of this class is their high 
ratio of m/m. Indeed high luminosity AGN which do not ful- 
fil the velocity width criteria of NLSls (i.e. <2000km/s) but 
share many of their other charac teristics are also th ought 
to have high m/m, e.g. PDS456 (iReeves et al.|[20o3 ). con- 
sistent with the model of lMcHardv et alH 20061 ). It is partly 
for this reason that NLSls have attracted much interest. 
In addition, we can potentially test the full extent of the 
relationships between observable parameters that have been 
established for the more general AGN population, such as re- 
lations between black hole mass a nd stellar bulge luminosity 
and velocity dispersio n (see e.g. iFerrarese fc Merrit3 12OO0I : 
iGebhardt et~ai]|2000l ). Also, evidence of some of the more 
dramatic features of AGN activity, in particular relativis- 
tic and non-relativistic outflows, arc sometimes prominent 
in NLSl spectra (e.g. Lcighlv 2001; Lcighly & Moore 20oJ), 
although it is unclear whether this is an indication of a real 
propensity toward outflows in NLSls, or whether other spec- 
tral features simply allow these features to be more easily 
measured in these AGN. This has encouraged their study as 
laboratories for some phenomena present in general in AGN, 
but exhibited in a more extreme form in the NLSls. 

The optical spectra of NLSl provide a wealth of in- 
formation. The large number of permitted and forbid- 
den lines give information on the kinematics and physical 
properties of the various emitting regions, as well as the 
shape of the (unobsc ured) ionising continuum. For example, 
iDietrich et al.l (|2005l ) compared the low ionisation (<100eV) 
emission lines of NLSl to show that the far-UV continuum 
of NLSl was not significantly different from that of BLSls, 
thereby indicating that the soft X-ray excess, seen in many 
NLSls, is unlikely to extend to significantly lower energies 
approaching lOOeV. However the true extent of this con- 
tinuum feature, which is likely to contribute a significant 
fraction of the bolometric luminosity in a large proportion 
of NLSls, is still unknown. 

With ionisation potentials greater than lOOeV, FHILs 
offer the possibility of determining the true shape of the 
soft X-ray excess. Previous studies of the optical FHILs 
in NLSls, most notably [Fe VIIJA6087 (ionisation energy: 
99eV), [Fe X]A6374 (235eV), [Fe XI]A7892 (262eV) and 
[Fe XIV]A5303 (361eV) noted that these lines were signif- 
icantly blueshifted and broadened, with respect to the low- 
ionisation forbidden lines, suggesting the possibility of out- 
flows between the traditi o nal broad- and narrow- line regions 
(e.g. lErkens et all Il997l : iPorguet et al] Il999l : iNagao et all 
I2OO0I ). However, despite the fact they often display broad 
wings and have prominent profile asymmetries, little effort 
has been directed towards modelling the FHILs using mul- 
tiple component fits. This could yield insights into the de- 
tailed kinematics of the highly ionised gas close to the AGN. 
The information gained from gas kinematics would help us 
to constrain the location of the FHIL emitting gas with re- 
spect to other emitting regions, and also help us obtain a 
better understanding of outflows from AGN. 

In this paper, we examine the optical emission line spec- 
tra of 10 nearby (z < 0.1) AGN at sufficiently high spectral 
resolution and signal to noise to permit us to model their 
FHILs using multiple components. Six of these AGN have 
previously been classified as NLSls. We also observed three 
BLSls to facilitate comparison between the Seyfert groups. 



In an attempt to determine whether the FHILs are affected 
by the obscuration of the putative dusty molecular torus, we 
also observed one Seyfert 2 galaxy, Mrk573, which has been 
shown to contain a hidden broad line region bv jNagao et al.l 
(|2004 ). The high quality of our spectra also enable accurate 
multiple component fits to both the forbidden and permitted 
emission lines, allowing precise measurement of the kinemat- 
ics of the broad and narrow emission line regions in NLSls 
and BLSls. This paper is structured as follows: in §2 we 
discuss the observations and spectral reduction techniques. 
We present our main results, including measurement of the 
FHILs, in §3 and in §4 we interpret our findings in terms 
of current AGN models. We conclude with a summary our 
principal results in the final section. 



2 OBSERVATIONS 

The ten AGN in our sample were observed using the ISIS 
long slit spectrograph on the William Herschel Telescope 
during the nights of the 17*'' and 18"" of October 2006. Ba- 
sic information on these objects is listed in Table [T] Both 
the red and blue arms of the ISIS instrument were used, 
with a 600 lines/mm grism in each case. The CCD detectors 
used were the EEV12 and the Marconi CCD in the blue and 
red arms respectively. In order to cover the wavelengths of 
all four principal optical FHILs ([Fe VII], [Fe X], [Fe XI], 
[Fe XIV]) using this instrument, it was necessary to observe 
each galaxy with two different grating angles in one of the 
arms. Based on pre vious studies of FHILs in NLSl galaxies 
(|Erkens et al.lll997l '). it was apparent that the [Fe XIV]A5303 
line was likely to be the weakest species of those observed. 
For this reason this line was observed twice using the blue 
arm, whilst the grating angle in the red arm was initially set 
for the measurement of the [Fe VII]A6087 and [Fe X]A6374 
species, followed by a grating angle change to allow mea- 
surement of the [Fe XI]A7892 line at a longer wavelength. 
This procedure thus resulted in a longer integration time 
for the weaker [Fe XIV] line than for the other iron species. 
The choice of beam-splitting dichroic was based on the ob- 
served wavelengths of the FHILs. Objects with z ^ 0.035 
were observed using the 5700A dichroic, while objects of 
z > 0.035 were observed using the 6IOOA dichroic. This 
choice ensured that the lines of most interest always lay in a 
wavelength region of high dichroic throughput. The instru- 
mental set-up provided almost complete spectral coverage 
between 4800A and 8900A. In order to minimise the stel- 
lar contribution from the host galaxy we used a slit width 
of 0.5", resulting in an instrumental spectral resolution of 
R = A/AA = c/Au ~ 3500 — 5500, depending on wave- 
length (equivalent to 85-55 km/s FWHM). Corresponding 
standard star and telluric dividing star frames were taken 
for each galaxy observation. Wavelength calibration frames 
were taken using a CuNe + CuAr arc lamp. 



2.1 CCD Image Reduction 

The raw CCD images first had the bias level subtracted, 
then they were flat fielded and trimmed using the IRAF 
set of image reduction software. A two dimensional disper- 
sion correction and wavelength calibration function was de- 
rived from the calibration arc frames. This was applied to 
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Galaxy 


Type 


z 


R.A. 


Dec 




Exposure Times (s) 


Ref. 








(J2000) 


Blue 


Red/Short 


Red /Long 




Ark 564 


NLSl 


0.0247 


22:42:39.3 


+29:43:31 


1080 


540 


540 


1 


1H1934-063 


NLSl 


0.0106 


19:37:33.0 


-06:13:05 


1800 


900 


900 


1 


IIZW136 


NLSl 


0.0633 


21:32:27.8 


+10:08:19 


2400 


1200 


1200 


2,3 


IZWl 


NLSl 


0.0611 


00:53:34.9 


+12:41:36 


2400 


900 


1200 


4 


Mrk335 


NLSl 


0.0261 


00:06:19.5 


+20:12:10 


1600 


600 


1000 


1,6 


Mrk573 


Sy2 


0.0172 


01:43:57.8 


+02:21:00 


1800 


900 


900 


5 


Mrk618 


NLSl 


0.0354 


04:36:22.2 


-10:22:34 


3000 


1500 


1500 


6 


NGC985 


Syl 


0.0432 


02:34:37.8 


-08:47:15 


1800 


1200 


1600 


7 


NGC7469 


Syl 


0.0163 


23:03:15.6 


+08:52:26 


900 


450 


450 


8 


VIIZW118 


Syl 


0.0803 


07:07:13.1 


+64:35:59 


4000 


2000 


2000 


9 



Table 1. AGN observed using ISIS on the WHT. Targe t co- ordinates taken from NA SA /IPAC Extra galactic Database. S pectral clas- 
sificat ions taken from: 1. iRodrfguez-Ardila et al.l ll2002f). 2. iBoroson fc Greeiil lll992l'). 3lConstantin fc Shields (200 31). 4. iRudv et al] 
1I2OOOI'). 5. NASA/IPAC Ex tragalactic Database, 6. lRvan et al.l | |2007^ . 7. lde Vaucouleurs fc de VaucouleursI 019751 ) 8. IRiffel et al.l 1I2OO6I ) 
9. iKunth fc Sargent! l ll979l) . Redshifts measured from the [SII] AA6717/32 doublet (assumed to represent the systemic velocity of the 
galaxy), with heliocentric correction. A typical error on these redshifts is ~45km/s 



each science frame, again using standard /i?ylF routines. The 
quality of the wavelength calibration was checked and cor- 
rected w here necessary us i ng th e sky emission lines iden- 
tified by lOsterbrock et al.l (|l996l ). To extract the ID spec- 
trum from the 2D image, an extraction routine was writ- 
ten using the IDL programming language. We initially col- 
lapsed the 2D image along the direction of dispersion in 
order to obtain an indication of the spatial profile of the 
object within the slit. This 'average' profile was then fit- 
ted with a single Gaussian. We then returned to the origi- 
nal 2D image and extracted the ID spectrum by summing 
the contribution of each pixel within 2a limits of the av- 
erage Gaussian profile, thus ensuring the same extraction 
window was used along the spectrum. The underlying sky 
background was determined by extrapolation of the region 
between Scr and 20a either side of the object spectra. These 
spectra were then fiux calibrated using the standard star ob- 
servations and the IRAF routines noao.onedspec. standard, 
noao.onedspec.sensfunc and noao.onedspec. calibrate. An at- 
tempt was then made to correct for telluric absorption by 
the division of a normalised, continuum removed (i.e. fiat) 
spectrum of a standard star. However, since this introduces 
large uncertainties in the spectra between the wavelengths 
affected by telluric absorption features, we have not at- 
tempted to obtain strengths for coronal lines within these 
regions. 

2.2 General Corrections 

Due to the procedure used to obtain the required wavelength 
coverage, there are three associated spectra for each galaxy; 
blue (4400-5600A), red/short (5700-7000A) and red/long 
(73OO-86OOA), which we have joined together to form a sin- 
gle, almost continuous spectrum. During the spectral joining 
process, it was noted that no significant corrections were 
needed to force the spectra covering the three wavelength 
ranges to align, thus giving us confidence in the flux calibra- 
tion procedure and the photometric accuracy. All our spec- 
tra have been corrected for redshift, which was calculated 
from the average observed wavelength of the [S II]A6717, 
6732 doublet, and geocentric velocities have been taken into 
account. All velocities published here are therefore relative 



to the [S II]A6717, 6732 doublet, which we assume is repre- 
sentative of the host galaxy. This is based on the fact that 
it has a low excitation potential and low critical density, 
meaning that it is likely to be formed away from the AGN. 

The spectra have been corrected for Galactic extinc- 
tion using the wavel ength depend ent absorption f unction of 
ICardeUi et ^ (|l989f ) and the Schl egel et alfdiggsl ) E(B-V) 
values obtained from NASA/IPAC Extragalactic Database. 
Due to the narrow slit widths used, we expect that the spec- 
tra suffer only minor contamination from the host galaxy. To 
test this assumption, we searched for the stellar absorption 
feature Mg b at A5170A, but only found evidence of it in 
Mrk573 (the only Seyfert 2 in our sample), in which it has 
an equivalent width of ~5A . Because we concentrate on 
the profiles and relative strengths of emission lines at wave- 
lengths which are largely unaffected by the stellar features, 
we have not corrected for this contamination. 

2.3 Continuum Fitting 

Since we wish to measure the total flux and profile of the 
emission lines, it was necessary to accurately fit and subtract 
the underlying continuum. For blends of weak lines, this 
can be problematic as it is difficult to distinguish between 
the continuum and extended line wings, thereby making the 
continuum fit somewhat subjective. In an attempt to more 
objectively determine the continuum levels, we modelled the 
whole spectrum with a gently inflecting fourth order poly- 
nomial fitted to regions of the spectrum that were perceived 
by eye to be free of emission features. Each fit was then in- 
spected at the wavelength of the lines of interest to check 
that it was a sensible fit to the local underlying continuum. 
As some of the weaker emission lines lie on the broad wings 
of the Hq profile, we first fit the Ha line and remove the 
contribution of the broad wings from the weak lines. 

2.4 Fe II Lines 

Permitted emission from Fe II is often detected in the optical 
spectra of AGN, and is particularly strong in many NLSls, 
in which multiplets of broad lines blend together to form a 
pseudo-continuum. As some of our lines of interest lie close 
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Figure 1. Flux calibrated, rcdshift corrected, Fell subtracted spectra. Flux has been normalised to the level at 5100Aand offset for 
clarity. The dotted lines either side of the plot show the level of the offset. For a clearer depiction of the lines discussed here, we refer 
the reader to Figs. |A31- El0l 



to strong Fe II features, we decided to rem ove them by using 
a met hod closely based on that outlined in lBoroson fc GreenI 
l|l992h . In their study of the strong Fe II hues in the spec- 
trum of IZWl, IVeron-Cettv et al'l (120041 ) show that this 
emission contains two velocity components. To account for 
this we used their published line lists as a template, and 
broadened the two sets of lines by convolving this template 
with two Gaussians with independent variables of veloc- 
ity wi dth and intensity. In their method, iBoroson fc GreenI 
l| 19921 ) produced a matrix of Fe II spectra, from which they 
selected the one that most closely matched the Fe II emis- 
sion in a particular object's spectrum. Increasing comput- 
ing processing speeds since their original work fortunately 
means that we can now produce the broadened Fe II spec- 
trum 'on the fly' using an IDL based GUI REMOVEFEII 
0. This enables much more control over the range of input 
variables, making the Fe II feature rem oval a more dynami c 
and controlled process. As was noted bv lLandt et"al] (|200'i1 ). 
it is difficult to determine the true widths of the Fe II lines as 
above a certain velocity width the lines blend together into 
a single broad feature. To counter this it was noted that 
by fitting a single line that was unique to either the nar- 
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row or broad component templates and that was not part 
of an Fe II blend, a reasonably independent estimate of the 
strength and width of the two templates could be deter- 
mined. A limitation to this method is that it is dependent 
on the assumption that all our AGN display the same Fe II 
emission lines, with the same relative strengths, as found in 
IZWl. 

In practice it was found that, in general, this procedure 
was able to remove a significant proportion of the narrow Fe 
II emission component, but it was difficult to fully account 
for the broad component of the emission lines. A qualita- 
tive assessment of the line profiles affected by Fe II emission 
is therefore possible, but the level of uncertainty does not 
allow reliable measurements in cases where the Fell blends 
are strong -between AA4100-4600 and AA5150 -5400 (in par- 
ticular the region of [Fe XIV]A5303). 

2.5 Line fitting 

All the emission lines discussed in detail in the subsequent 
sections were modelled by fitting multiple Gaussian pro- 
files, which we argue is more statistically robust than the 
Lorentzian component fitting metho d used in some previ- 
ous NLSl studies (see o ur appendix; lGon9alves et al.lll999l . 
IVeron-Cettv et al.llioOll '). Blends (excluding the Fe II lines 
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Figure 1 - continued 



discussed above) were dealt with by simultaneously fitting 
the blended lines. This method worked well for the majority 
of blended lines. In regions of the spectra where broad lines 
are blended together with narrow lines, (in particular the 
Ha and [NII]AA6548,6584 lines), once we had fitted these 
lines we compared the widths of the blended narrow lines 
with those of the unblended narrow [OIII]AA4959,5007 lines 
as a check that the fit that been performed correctly. 

The 90% confidence errors associated with the param- 
eters derived for each fitted Gaussian (i.e. intensity, width 
and central wavelength) were determined by deviating each 
of the parameters from their optimum values until a change 
in equal to the published 90% confidence level for the 
appropriate number of free fitting parameters was obtained. 
In addition to these random errors, we found that uncer- 
tainties in the continuum model leads to significant system- 
atic errors for the line parameters (especially the broadest 
components). To estimate this, we rescaled the continuum 
by ±lcr (adopting RMS as la) and refit the line models, 
thereby quantifying the impact of continuum levels that are 
set too high or too low. 

The intensities, widths and shifts of the multiple Gaus- 
sian fits for the lines discussed in this paper, along with their 
associated random and systematic errors, are listed in Tables 
[2]-|4l which we discuss in turn below. We note that, in Ta- 
ble |3] we have not corrected for instrumental broadening as 



none of the components of the lines of interest had widths 
close to the instrumental broadening determined from the 
sky lines (~2A). Therefore, when considered in quadrature, 
instrumental broadening would have an insignificant effect 
on the measured line widths. 



3 EMISSION LINE RESULTS 
3.1 The Broad Emission Lines 

Due to the clear distinction between the spectrum of the 
only Seyfert 2 in our sample, Mrk573, and all the other 
galaxies we observed (Seyfert Is, with a range of broad emis- 
sion line widths), we will discuss this galaxy in a separate 
section. 



3.1.1 Balmer Lines 

As expected, all our spectra are dominated by strong Ha 
and H/3 lines, which we show, together with their Gaussian 
fits, in Figs IA3l and IA4| respectively. We have been able to 
deconvolve these lines into their narrow and broad compo- 
nents for all the observed Seyfert 1 nuclei. For six of the 
nine Seyfert Is we observed the Balmer lines are best fitted 
with three components. The three galaxies that do not fit 
this trend are NGC7469, VIIZW118 and NGC985, in which 
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Table 2. Modelled line flux, as a percentage of the H/3 line. Where no value is listed, there was no significant indication of the presence 
of the line or component. The ± errors refer to the 90% confidence interval of the measurement, the (%) errors refer to the systematic 
errors associated with the estimation of the continuum fiux levels. The H/3 line flux has units of 10~^^ ergs/s/cm^ 



Line 


Ark564 


IZWl 


1H1934-063 


Mrk335 


NGC7469 


np Flux 


2.76 


2.62 


6.11 


6.70 


5.43 



HaA6563 
Broad 

Intermediate 
Narrow 
H^A4861 
Broad 

Intermediate 

Narrow 
[O III]A5007 

Broad 

Narrow 
He IA5876 

Broad 

Narrow 
He IIA4686 

Broad 

Narrow 
[Fe VII]A6087 

Broad 

Narrow 
[Fe X]A6374 

Broad 

Narrow 
[Fe XI]A7892 

Broad 

Narrow 
[Fe XIV]A5303 

Broad 



454.2 ± 5.0 (5%) 
184.9 ± 2.5 (<2.5%) 

169.1 ± 2.1 (10%) 

100.2 ± 1.3 (5%) 
100.0 ± 1.4 (10%) 

46.80 ± 0.92 (5%) 
36.36 ± 0.58 (10%) 
16.84 ± 0.33 (20%) 
88.50 ± 0.96 (5%) 

44.02 ± 0.54 (<2.5%) 
44.48 ± 0.51 (5%) 
17.0 ± 1.1 (45%) 

12.81 ± 0.99 (60%) 
4.15 ± 0.37 (<2.5%) 
25.54 ± 0.64 (15%) 
18.97 ± 0.57 (10%) 
6.57 ± 0.25 (25%) 
3.95 ± 0.29 (25%) 

3.95 ± 0.29 (25%) 

9.36 ± 0.66 (25%) 
8.03 ± 0.56 (25%) 
1.32 ± 0.33 (25%) 
8.83 ± 0.57 (30%) 

5.96 ± 0.49 (30%) 
2.88 ± 0.30 (35%) 
1.81 ± 0.33 (35%) 
1.81 ± 0.33 (35%) 



453.8 ± 3.3 (15%) 
228.0 ± 1.9 (10%) 
175.5 ± 1.4 (5%) 
50.33 ± 0.59 (50%) 
100.00 ± 0.92 (10%) 
61.66 ± 0.69 (15%) 
28.18 ± 0.33 (5%) 
10.16 ± 0.19 (10%) 
65.10 ± 0.61 (15%) 
48.66 ± 0.49 (15%) 
16.45 ± 0.26 (10%) 



320.4 ± 1.5 (5%) 
76.06 ± 0.78 (<2.5%) 

198.03 ± 0.86 (5%) 
46.27 ± 0.38 (<2.5%) 
100.00 ± 0.48 (15%) 
43.10 ± 0.30 (10%) 
50.75 ± 0.25 (20%) 
6.15 ± 0.10 (15%) 
38.87 ± 0.21 (<2.5%) 
28.19 ± 0.17 (<2.5%) 
10.679 ± 0.098 (<2.5%) 
12.22 ± 0.34 (35%) 
8.59 ± 0.30 (40%) 
3.63 ± 0.17 (25%) 
22.52 ± 0.24 (15%) 
15.26 ± 0.21 (5%) 

7.25 ± 0.11 (35%) 
3.44 ± 0.14 (30%) 
3.44 ± 0.14 (30%) 

6.26 ± 0.21 (25%) 
6.26 ± 0.21 (25%) 

7.80 ± 0.23 (35%) 

7.80 ± 0.23 (35%) 

1.81 ± 0.11 (60%) 
1.81 ± 0.11 (60%) 



350.1 ± 2.4 (<2.5%) 
155.0 ± 1.4 (<2.5%) 
139.4 ± 1.1 (<2.5%) 
55.77 ± 0.66 (5%) 
100.00 ± 0.74 (5%) 
60.16 ± 0.52 (5%) 
28.92 ± 0.30 (5%) 
10.92 ± 0.20 (5%) 
29.32 ± 0.25 (5%) 
18.19 ± 0.20 (5%) 
11.13 ± 0.12 (5%) 
23.27 ± 0.81 (30%) 
16.45 ± 0.74 (35%) 
6.82 ± 0.32 (15%) 
12.18 ± 0.29 (10%) 
8.07 ± 0.25 (5%) 

4.11 ± 0.15 (20%) 

3.12 ± 0.26 (65%) 
3.12 ± 0.26 (65%) 

5.68 ± 0.59 (55%) 
4.80 ± 0.54 (60%) 
0.88 ± 0.24 (15%) 
2.24 ± 0.17 (60%) 
2.24 ± 0.17 (60%) 



487.5 ± 3.2 (5%) 
424.1 ± 2.8 (5%) 

63.41 ± 0.64 (<2.5%) 
100.00 ± 0.83 (5%) 
85.76 ± 0.74 (5%) 

14.24 ± 0.24 (<2.5%) 
111.13 ± 0.77 (<2.5%) 

46.32 ± 0.44 (5%) 
64.80 ± 0.45 (<2.5%) 

23.33 ± 0.89 (35%) 
18.05 ± 0.80 (40%) 
5.27 ± 0.38 (10%) 
10.71 ± 0.39 (40%) 
8.83 ± 0.36 (45%) 
1.88 ± 0.15 (5%) 
5.35 ± 0.66 (40%) 
2.79 ± 0.55 (50%) 
2.56 ± 0.35 (30%) 
7.29 ± 0.68 (70%) 
7.29 ± 0.68 (70%) 



Line 


Mrk618 


IIZW136 


VIIZW118 


NGC985 


Mrk573 


H/3 Flux 


1.84 


4.64 


1.47 


3.63 


0.400 



H<:«A6563 
Broad 

Intermediate 
Narrow 
H^A4861 
Broad 

Intermediate 

Narrow 
[O III]A5007 

Broad 

Narrow 
He IA5876 

Broad 

Narrow 
He IIA4686 

Broad 

Narrow 
[Fe VII]A6087 

Broad 

Narrow 
[Fe X]A6374 

Broad 

Narrow 
[Fe XI]A7892 

Broad 

Narrow 
[Fe XIV]A5303 

Broad 



331.1 ± 3.4 (5%) 
184.1 ± 2.2 (5%) 
136.5 ± 1.6 (10%) 
10.54 ± 0.42 (<2.5%) 
100.0 ± 1.2 (20%) 
76.37 ± 0.95 (10%) 
22.03 ± 0.48 (45%) 
1.61 ± 0.24 (40%) 
80.54 ± 0.78 (5%) 
65.06 ± 0.64 (5%) 
15.48 ± 0.23 (<2.5%) 
26.0 ± 1.1 (35%) 
18.8 ± 1.0 (40%) 
7.22 ± 0.53 (25%) 
18.83 ± 0.63 (20%) 
13.51 ± 0.57 (25%) 
5.32 ± 0.26 (5%) 
4.29 ± 0.30 (35%) 
4.29 ± 0.30 (35%) 

6.60 ± 0.45 (30%) 
6.60 ± 0.45 (30%) 



383.0 ± 2.3 (5%) 
116.2 ± 1.00 (<2.5%) 

209.7 ± 1.3 (5%) 
57.15 ± 0.47 (<2.5%) 
100.00 ± 0.76 (5%) 
34.50 ± 0.46 (5%) 
55.84 ± 0.42 (10%) 
9.66 ± 0.17 (<2.5%) 
23.94 ± 0.26 (5%) 
12.97 ± 0.21 (5%) 
10.97 ± 0.12 (5%) 



5.00 ± 0.24 (45%) 
3.00 ± 0.21 (65%) 
2.00 ± 0.11 (20%) 
4.27 ± 0.26 (40%) 
4.27 ± 0.26 (40%) 

6.47 ± 0.32 (45%) 
6.47 ± 0.32 (45%) 



412.8 ± 3.4 (10%) 
122.0 ± 1.6 (10%) 

290.8 ± 2.3 (10%) 
100.0 ± 1.0 (5%) 
41.11 ± 0.66 (10%) 

58.89 ± 0.60 (5%) 
20.54 ± 0.37 (10%) 
11.08 ± 0.28 (10%) 
9.46 ± 0.22 (5%) 
25.4 ± 1.1 (20%) 
17.48 ± 0.94 (25%) 
7.94 ± 0.51 (5%) 
5.28 ± 0.39 (35%) 
3.16 ± 0.33 (55%) 
2.12 ± 0.21 (10%) 
2.81 ± 0.31 (40%) 
2.81 ± 0.31 (40%) 

5.57 ± 0.59 (80%) 
5.57 ± 0.59 (80%) 



527.5 ± 1.9 (<2.5%) 
501.4 ± 1.8 (<2.5%) 

26.09 ± 0.32 (5%) 
100.00 ± 0.46 (10%) 
95.07 ± 0.43 (10%) 

4.926 ± 0.095 (10%) 
60.36 ± 0.27 (<2.5%) 
25.85 ± 0.18 (<2.5%) 
34.51 ± 0.15 (<2.5%) 
15.68 ± 0.62 (40%) 
14.77 ± 0.59 (40%) 
0.91 ± 0.18 (15%) 
1.145 ± 0.050 (25%) 
1.145 ± 0.050 (25%) 

2.80 ± 0.15 (50%) 
2.80 ± 0.15 (50%) 



382. ± 10.0 (<2.5%) 

242.8 ± 6.8 (<2.5%) 

138.9 ± 4.0 (<2.5%) 
100.0 ± 3.3 (10%) 

66.0 ± 2.5 (5%) 
34.0 ± 1.4 (15%) 
1216. ± 28. (<2.5%) 

785. i 18. (<2.5%) 
431. ± 10. (<2.5%) 
9.1 ± 1.4 (45%) 
9.1 ± 1.4 (45%) 

35.6 ± 1.9 (20%) 
22.0 ± 1.5 (25%) 

13.7 ± 1.0 (10%) 
28.0 ± 1.9 (20%) 

28.0 ± 1.9 (20%) 

15.4 ± 2.1 (30%) 
15.4 ± 2.1 (30%) 

11.1 ± 1.4 (80%) 
11.1 ± 1.4 (80%) 
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Table 3. Line widths in km/s. As above, the random errors are the 90% confidence interval and the systematic uncertainties arising 
from the continuum model are given in parenthesis, measurement. 



Line 


Ark564 




IZWl 


1H1934-063 


Mrk335 


NGC7469 


HaA6563 
















Broad 


2578. ± 24. (10%) 


3247. 


± 


16. (20%) 


2837. ± 26. (10%) 


3920. ± 27. (5%) 


2066.1 ± 6.4 (<2.5%) 


Intermediate 


963.1 ± 7.3 (5%) 


918.3 


± 


4.1 (25%) 


1271.7 ± 3.7 (<2.5%) 


1333.4 ± 7.3 (<2.5%) 


- 


Narrow 


320.9 ± 3.2 (<2.5%) 


328.2 


± 


3.2 (20%) 


411.4 ± 3.2 (<2.5%) 


707.1 ± 8.2 (<2.5%) 


230.4 ± 2.3 (<2.5%) 


H^A4861 
















Broad 


3175. ± 60. (25%) 


4049. 


± 


43. (15%) 


2624. ± 18. (25%) 


3838. ± 29. (<2.5%) 


2046. ± 15. (5%) 


Intermediate 


897. it 12. (15%) 


967. 


± 


11. (5%) 


1059.7 ± 4.3 (10%) 


1292. ± 12. (<2.5%) 


- 


Narrow 


290.1 ± 5.6 (10%) 


388.8 


1 ± 


7.4 (5%) 


317.8 ± 6.2 (10%) 


685. ± 14. (<2.5%) 


329.6 ± 6.2 (<2.5%) 


[O III]A5007 
















Broad 


336.1 ± 2.4 (5%) 


2044. 


± 


20. (5%) 


538.0 ± 3.0 (<2.5%) 


893. ± 11. (5%) 


744.2 ± 6.6 (10%) 


Narrow 


161.2 ± 1.2 (<2.5%) 


759. 


± 


13. (5%) 


191.1 ± 1.8 (<2.5%) 


246.3 ± 2.4 (5%) 


284.6 ± 1.2 (<2.5%) 


He IA5876 
















Broad 


3250 ± 510 (25%) 






- 


1852. ± 72. (45%) 


6970 ± 380 (30%) 


2910 ± 160 (25%) 


Narrow 


475. ± 46. (<2.5%) 








576. ± 31. (15%) 


1336. ± 76. (10%) 


665. ± 57. (5%) 


He IIA4686 
















Broad 


2305. ± 76. (20%) 








2535. ± 37. (20%) 


2666. ± 93. (5%) 


1700 ± 110 (25%) 


Narrow 


474. ± 22. (20%) 








686. ± 12. (20%) 


957. ± 40. (10%) 


290. ± 25. (5%) 


[Fe VII]A6087 
















Broad 


446. ± 37. (15%) 








572. ± 27. (20%) 


1380 ± 140 (45%) 


1270 ± 360 (45%) 


Narrow 














518. ± 85. (20%) 


[Fe X]A6374 
















Broad 


749. ± 62. (25%) 








1127. ± 42. (15%) 


2680 ± 420 (35%) 


2560 ± 300 (35%) 


Narrow 


268. ± 73. (10%) 










510 ± 170 (10%) 




[Fe XI]A7892 
















Broad 


1140 ±110 (40%) 








1457. ± 48. (20%) 


1330 ± 120 (35%) 




Narrow 


428. ±51. (20%) 














[Fe XIV]A5303 
















Broad 


540 ±110 (30%) 








623. ± 55. (35%) 







Line Mrk618 IIZW136 VIIZW118 NGC985 Mrk573 
HaA6563 

Broad 3221. ± 26. (5%) 4474. ± 30. (5%) 5242. ± 59. (20%) 4987.5 ± 7.8 (<2.5%) 

Intermediate 1681. ± 15. (<2.5%) 2093.1 ± 5.9 (<2.5%) - - 533.9 ± 9.1 (<2.5%) 

Narrow 191.1 ± 8.2 (<2.5%) 675.1 ± 4.6 (<2.5%) 2503.1 ± 8.7 (<2.5%) 519.3 ± 8.7 (5%) 180.1 ± 3.7 (<2.5%) 
H^A4861 

Broad 3343. ± 35. (10%) 5130. ± 68. (20%) 4543. ± 70. (20%) 5003. ± 18. (5%) 

Intermediate 1360. ± 33. (30%) 2353. ± 14. (5%) - - 546. ± 19. (20%) 

Narrow 240. ± 64. (5%) 649. ± 12. (<2.5%) 2293. ± 18. (<2.5%) 467. ± 12. (10%) 185.8 ± 7.4 (5%) 
[O III]A5007 

Broad 810.1 ± 4.8 (5%) 1379. ± 26. (15%) 772. ± 22. (10%) 803.5 ± 5.4 (5%) 523.1 ± 2.4 (<2.5%) 

Narrow 267.2 ± 3.6 (<2.5%) 339.1 ± 3.6 (5%) 442. ± 11. (<2.5%) 309.8 ± 1.2 (<2.5%) 169.6 ± 1.2 (<2.5%) 
He IA5876 

Broad 5500 ± 340 (35%) - 5470 ± 360 (15%) 5150 ± 280 (25%) 254. ± 45. (30%) 

Narrow 1600 ± 130 (10%) - 1810 ± 140 (<2.5%) 520 ± 130 (10%) 
He IIA4686 

Broad 4110 ± 190 (10%) 2990 ± 270 (45%) 2070 ± 260 (25%) 423. ± 22. (15%) 430. ± 30. (20%) 

Narrow 920. ± 50. (5%) 1061. ± 68. (5%) 940 ± 100 (10%) - 198. ± 17. (5%) 
[Fe VII]A6087 

Broad 775. ± 63. (15%) 2020 ± 190 (20%) 1060 ± 140 (25%) 980. ± 65. (35%) 378. ± 30. (15%) 

Narrow _____ 
[Fe X]A6374 

Broad 1290 ± 130 (15%) 2740 ± 180 (25%) 2980 ± 400 (45%) - 387. ± 58. (20%) 

Narrow _____ 
[Fe XI]A7892 

Broad _ _ _ _ 330. ± 43. (45%) 

Narrow _____ 
[Fc XIV]A5303 

Broad _____ 
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Table 4. Line shift from systemic velocity shift (in km/s), as determined from the [S II]A6718 doublet. We use the convention that 
negative numbers indicate a net redshift. 



Line 



Ark564 



IZWl 



1H1934-063 



Mrk335 



NGC7469 



H„A6563 

Broad 

Intermediate 
Narrow 
H^A4861 
Broad 

Intermediate 
Narrow 
[O III]A5007 

Broad 

Narrow 
He IA5876 

Broad 

Narrow 
He IIA4686 

Broad 

Narrow 
[Fe VII]A6087 

Broad 

Narrow 
[Fe X]A6374 

Broad 

Narrow 
[Fe XI]A7892 

Broad 

Narrow 
[Fe XIV]A5303 

Broad 



-71. ± 13. 
-32.4 ± 4.1 
-7.3 ± 1.8 

-35. ± 31. 
-7.3 ± 6.8 
-10.4 ± 3.1 

13.6 ± 1.8 
-42.12 ± 0.60 

-220 ± 180 
-8. ± 25. 

-23. ± 38. 
55. ± 11. 

73. ± 19. 



256. ± 31. 
86. ± 41. 

424. ± 55. 
170. ± 26. 

425. ± 59. 



234.4 ± 8.2 
99.6 ± 2.3 
-19.2 ± 1.8 

461. ± 21. 

252.1 ± 5.6 
72.0 ± 3.7 

1375.8 ± 9.6 

537.2 ± 6.6 



-108. ± 14. 
10.5 ± 1.8 
-0.5 ± 1.8 

-76.3 ± 9.9 
0.8 ± 2.5 
-63.9 ± 3.1 

105.9 ± 1.8 
-37.9 ± 1.2 

-369. ± 38. 
52. ± 16. 

30. ± 19. 
133.5 ± 6.4 

179. ± 14. 



118. ± 14. 
298.3 ± 0.00 
-8.4 ± 4.6 

90. ± 12. 
257.5 ± 6.2 
37.9 ± 6.2 

149.0 ± 5.4 
-15.1 ± 1.2 

-230 ± 190 
141. ± 37. 

-428. ± 45. 
161. ± 19. 

359. ± 69. 



510. ± 19. 450 ± 190 

228. ± 83. 

396. ± 23. 473. ± 60. 



720. ± 24. 



-249.2 ± 0.00 

-81.5 ± 0.00 

-312.7 ± 6.2 

-91.9 ± 0.00 

215.5 ± 3.6 
-20.44 ± 0.60 

-754. ± 79. 
-81. ± 29. 

59. ± 45. 
35. ± 13. 

520 ± 150 
99. ± 39. 

440 ± 150 



Line 



Mrk618 



IIZW136 



VIIZW118 



NGC985 



Mrk573 



HaA6563 
Broad 

Intermediate 
Narrow 
H^A4861 
Broad 

Intermediate 

Narrow 
[O III]A5007 

Broad 

Narrow 
He IA5876 

Broad 

Narrow 
He IIA4686 

Broad 

Narrow 
[Fe VII]A6087 

Broad 

Narrow 
[Fe X]A6374 

Broad 

Narrow 
[Fe XI]A7892 

Broad 

Narrow 
[Fe XIV]A5303 

Broad 



20. ± 14. 
-5.1 ± 8.2 
-96.9 ± 4.6 

125. ± 22. 

85. ± 17. 
-84. ± 54. 

235.2 ± 2.4 
-49.3 ± 2.4 

-410 ± 170 
202. ± 68. 

-651. ± 97. 
261. ± 26. 

258. ± 30. 



387. ± 52. 



32. ± 16. 
-29.9 ± 3.2 
177.1 ± 2.3 

525. ± 35. 

44.1 ± 7.4 
189.0 ± 6.8 

222. ± 12. 
-0.3 ± 2.4 



-870 ± 130 
426. ± 36. 

691. ± 74. 



640. ± 80. 



-553. ± 32. 
-186.6 ± 5.0 
-704. ± 38. 

-126.5 ± 9.9 

106. ± 11. 
-89.3 ± 6.0 

-1520 ± 170 
-146. ± 68. 

-920 ± 130 
133. ± 57. 

279. ± 69. 
480 ± 180 



-88.9 ± 4.6 
-20.4 ±4.1 
-195.0 ± 9.3 

-16.7 ± 5.6 

104.1 ± 0.00 
9.4 ± 0.00 

-260 ± 130 
-27. ± 62. 

21.2 ± 6.4 
81. ± 30. 



-137.5 ± 4.6 
-57.6 ± 0.00 



-135.1 ± 6.2 
-100.0 ± 0.00 

-131.2 ± 1.2 
-69.51 ± 0.60 

-61. ± 23. 



-109. ± 17. 
-61.4 ± 9.0 

-73. ± 15. 



-78. ± 28. 



-64. ± 23. 
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a third component does not provide a significant improve- 
ment to the fit. In general, we find equivalent velocity width 
components (to within errors) in both the Ha and H/3 lines, 
giving us further confidence in the reality of each of the com- 
ponents. The only case of which this is not strictly true is the 
intermediate component of 1H1934-063, although here the 
discrepancy is <20%- typical of the size of the systematic er- 
rors on Balmer line components in other galaxies considered 
here. 

To assist in the discussion we have arranged our results 
in order of increasing Ha intermediate component velocity 
width, which we find most closely matches the crude FWHM 
of the whole line. Fitting the Balmer lines with Gaussians is 
complicated by the presence of a number of blended narrow 
and broad lines emitted by other species. In particular, the 
Ha line is blended with two often quite strong, narrow [N II] 
lines at A6548 and A6584. Fortunately, as they are narrow 
these lines are easily fitted and removed when fitting the Ha 
line. More difficult to account for is the broad He I A4922 
line that for ms the 'red shelf o f the H/3 line, and is discussed 
in detail bv lVeron et al.l ()2002l ). If, during the initial fitting, 
there was evidence of a broad, redshifted component to the 
H/3 line, we assumed that it was the result of this Helium 
line. We then refit the blend, adding an additional Gaussian 
with its central wavelength fixed to 4922A(rest frame). This 
component is then excluded from further analysis of the H/3 
line. A similar procedure was used to remove the less promi- 
nent red shelf on the Ha line, caused by the presence of He 
I A6678 emission. 

3.1.2 Balmer Ratios in the Velocity Components 

It is possible that emission from the BELR and NELR expe- 
rience diflerent amounts of reddening depending on the loca- 
tion of the dust within the emitting regions. It is also possi- 
ble that the physical conditions of opacity (self-absorption) 
and coUisional effects will lead to differences in the Balmer 
ratios across the various velocity components of the Balmer 
lines. To test this possibility we have measured the Ha/H/3 
ratios for both emission components (see Table [SJ. We find 
that for all our targets the degree of reddening calculated 
from the Balmer decrement increases with decreasing com- 
ponent velocity widths (to within systematic errors), which 
could imply a higher concentration of dust in the NLR 
(narrow line region), as has previously been suggested by 
iNetzer fc Laoil (|l993l '). 

We note that in 1H1934-063 (and possibly Mrk618) 
the broad component has an Ha/H/3 ratio that is less than 
the value predicted if there were no absorption whatsoever. 
Therefore it is clear that effects other than reddening are 
at play. However, we should be aware that the intensity 
of the broadest component is the most uncertain of all the 
fitted components as it is mainly determined by the weak 
wings on either side of the line centre. These wings are of- 
ten blended, especially in the case of H/3, and in addition 
are highly dependent on the fitting of the continuum level. 
In contrast, the intermediate and narrow components are 
better constrained by the profile of the line where the line 
flux is strongest. 

Our work confirms that the overall profiles of Ha and 
H/3 are very often not the same. This implies differences in 
one or all the parameters of reddening, density and tem- 



Table 5. Ho/H/J ratios. As we discuss in tlie text tlie systematic 
errors associated witli the continuum flux level dominate over 
the random errors, so here we only give the systematic errors (in 
parentheses as percentages). Many of the galaxies have different 
Balmer ratios for each velocity component, suggesting different 
amounts of absorption and / or physical conditions in each region. 



Galaxy 


Total 


Broad 


Interm. 


Narrow 


Ark564 


4.54 (15) 


3.95 (5) 


4.65 (20) 


5.95 (25) 


IZWl 


4.54 (25) 


3.70 (25) 


6.23 (10) 


4.95 (60) 


1H1934-063 


3.20 (20) 


1.76 (10) 


3.90 (25) 


7.52 (15) 


Mrk335 


3.50 (5) 


2.58 (5) 


4.82 (5) 


5.11 (10) 


NGC7469 


4.88 (10) 


4.95 (10) 




4.45 (<2.5) 


Mrk618 


3.31 (25) 


2.41 (15) 


6.20 (70) 


6.57 (50) 


IIZW136 


3.83 (10) 


3.37 (5) 


3.75 (15) 


5.92 (<2.5) 


VIIZW118 


4.13 (15) 


2.97 (20) 




4.94 (15) 


NGC985 


5.28 (10) 


5.27 (10) 




5.30 (15) 


Mrk573 


3.82 (10) 




3.68 (5) 


4.09 (15) 



perature, although this cannot be fully invest igated without 
coverage of other B almer lines (H7, H(5; e.g IPopovid l2003l . 
iLa Mura et"aDl2007l l. 

3.1.3 Helium Emission 

Our spectra also cover a number He lines, of which He IA5876 
and He II A4686 are two of the strongest. Plots of these lines, 
including their Gaussian fits, are shown in Figs lA5l and lA6l 
respectively. We have been able to measure both these emis- 
sion lines in seven of the nine observed Syl galaxies (the red- 
shifts of IZWl and IIZW136 place their He IA5876 between 
the ranges of the Blue and Red/Short spectral settings). Un- 
fortunately, He IIA4686 suffers from strong blending with 
the Fe II lines and the H/3 line just blueward and redward of 
it, respectively. Despite our attempts to remove the blended 
Fe II emission there is still some remnants of this procedure 
(seen either in emission, or in absorption caused by overcom- 
pensation). There also appears to be an unidentified broad 
line present in the blue wing for the majority of the ob- 
served He IIA4686 that does not seem to change after the 
Fe II removal. Where this unidentified line is blended with 
the He IIA4686 line, we have fit it with a single Gaussian 
and have not included it in the interpretation of the He II 
line intensity and width. Aware of these complications, we 
fit the He II line with multiple Gaussians where the addition 
of a second component provides a statistically better fit. 

When it is present the broad component is generally 
redshifted with respect to the narrow component of the line, 
however, this may be an artefact of line asymmetries caused 
by the Fe II removal. For the galaxies in which He IIA4686 
is detected its strength relative to H/3 varies considerably 
between galaxies, ranging from almost 25% in Ark564, to 
less than 2% in NGC985. 

More accurately measured is the He IA5876 line. It is 
well detected in all our spectra that cover the region, and is 
easily deconvolved into its constituent velocity components. 
In 4/9 of the Syl-type AGN the broadest of these compo- 
nents is considerably broader than the broadest H/3 velocity 
component, although the level of the continuum fit has a 
considerable effect on the width of the fitted broad compo- 
nent (see 
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3.2 Narrow Emission Lines 

The narrow emission line region lias been well studied in 
the past, and is therefore much better understood than the 
other two non-spatially resolved regions discussed in this 
paper. Because of this we do not provide a full analysis of 
the narrow emission lines. Instead we will simply mention 
the main points of the most prominent of these narrow lines. 

3.2.1 [O III] XX4959,5007 

All our spectra with the possible exception of IZWl, dis- 
play strong [O III] AA4959,5007 lines, whose profiles all show 
asymmetries or have low intensity broad wings that can- 
not be modelled by a single Gaussian. These asymmetries 
tend toward a blue wing, which when fitted using two Guas- 
sians, suggests a narrow (200km/s<FWHM<400km/s) com- 
ponent centred on the systemic wavelength, plus a broader 
component (350km/s<FWHM<1000km/s) that is generally 
blueshifted with r espect to the nar r ow co mponent. This fact 
was first noted bv lHeckman et al] (ligsil ). in which they at- 
tribute the effect to outflows with preferentially absorbed 
redshifted components, supported by the tight correlation 
between line asymmetry and Ha/II/3 ratio. 

3.2.2 [S Ilj \6717, 6732 doublet 

As this feature is commonly used to determine the densities 
of the NLR of AGN, we chose to fit these lines with separate 
Gaussians, using multiple components if there was evidence 
of broad wings or inflections. Comparing the [S II] 6716/6731 
line ratio with published density relationships revealed no 
correlation between electron density and AGN type or width 
of H/3. This finding is consiste nt with the large scatter in the 
[SII] ratio vs. line width plot of lXu et al.l l|2007l ) at low values 
of H/3 FWHM. 

3.2.3 Balmer line narrow component 

As a test to determine whether the narrow component of the 
Balmer emission is arising from the same region of the AGN 
as the narrow forbidden line emission, we compare in detail 
the widths of the narrow components of the Balmer and [O 
III] lines. Table [3] shows that 7/9 of the observed Seyfert Is 
have narrow Balmer components that either lie in between 
or very close to the velocity widths of the narrow and broad 
components of the [O III] component widths. The two galax- 
ies in which this is not the case are IZWl and VIIZW118. 
As mentioned in the previous section, IZWl has very broad, 
highly blueshifted [QUI] lines, so it is not surprising that the 
widths do not match in this case, and there is no evidence 
of a narrow component in the Balmer lines of VIIZW118. 

3.3 Forbidden High lonisation (Coronal) Lines 
(FHILs) 

We detect the [Fe VII] A6087 emission line in every galaxy. 
The line profile in NGC7469 is best described by a two- 
Gaussian model. Additionally, Ark564 and Mrk335 both ap- 
pear to have a slight blue wing, but the improved fit of a 2 
component model is not statistically significant. 



We detect the [Fe X] emission line in seven of our spec- 
tra. Our measurements are complicated by blending with the 
[O I]A6363 line. An estimate of the contribution from the [O 
I] line can be made by using the nearby [O I]A6300 line as 
a calibration. From atomic physics these two lines should 
have exactly the same profile, and the intensities must have 
a ratio of 1:3. Using this information we are able to remove 
the contribution of [O I] from the [Fe X] blend. As with the 
He I and He II lines, measuring the broad component of the 
[Fe X] line is complicated by uncertainties in setting the lo- 
cal continuum levels and blends of multiple lines that form a 
pseudo-continuum. Despite these difficulties, we have mea- 
sured the [Fe X] lines using Gaussian fits, which are shown 
in Fig. IA8I Our analysis include a 2""* Gaussian component 
when statistically significant, although the majority of the 
spectra only require a single component. All the [Fe X] lines 
we measure contain broad components that are broader than 
the [Fe VII] broad components. The two galaxies that show 
the narrowest Balmer lines also have the narrowest [Fe X] 
lines of the sample. All the [Fe X] lines are blueshifted in the 
range between ~150 and ~700km/s, although there does not 
appear to be any correlation between blueshift and velocity 
width or relative intensity of the line between galaxies. 

Unfortunately, about half of our spectra have rather 
poor S/N ratios in the region containing the [Fe XI] line. 
This is partly attributed to telluric absorption (as discussed 
in i]2.ip . and partly to interference fringes on the CCD chip 
around these wavelengths. Due to this, we only have mea- 
surements of this line for three galaxies, Ark564, 1H1934-063 
and Mrk335, as shown in Fig. EE] The widths of the [Fe XI] 
lines reveals that they are similar to the [Fe X] line velocity 
widths in Ark564 and 1H1934-063. Mrk335 seems to have a 
wider broad component in [Fe X], but this value is uncertain 
because of poor continuum determination. 

The highest ionisation line species we observe, [Fe 
XIV]A5303, lies within a region that is heavily affected by 
the presence of Fe II emission. Despite our attempts at mit- 
igating the effects of Fe II emission, we have still had great 
difficulty in isolating the [Fe XIV]A5303 fine. The [Fe XIV] 
line is also strongly blended with the [Ca V]A5309 line. 
Although we cannot measure the [Fe XIV] line accurately 
enough to give robust quantitative results, we are confident 
of its presence in two of our sample, Ark564 and 1H1934-063 
(Ax'^ >500 and >1000, respectively, for 1 added Gaussian- 
corresponding to a 2>10(J detection in both targets). In both 
cases we are able to separate the feature into 2 Gaussian 
components (Ax^=52 and =85 between 1- and 2- compo- 
nent fit, again, corresponding to a >10a detection of two 
separate lines in both targets), which we interpret as [Fe 
XIV] and [Ca V]. We take further confidence in this from 
the fact that the velocity shift of the [Fe XIV] line is simi- 
lar to that of the [Fe X] and [Fe XI] lines. Differentiation of 
the [Fe XIV] and [Ca V] lines is helped by the fact that the 
former is more strongly blueshifted thereby increasing the 
observed separation from the intrinsic OA to 15-20A. That 
the ionisation potential of [Ca V] is significantly less than [Fe 
XIV] (67eV and 361eV, respectively) provides further evi- 
dence that blueshift increases with the ionisation potential 
of the emitting species. 
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3.4 Mrk 573 

Mrk573 is a relatively well studied nearby {z — 0.0170) 
Seyfert 2 galaxy, which has been shown t o contain a hidden 
BLR that is detectable in polarised light (|Nagao et al.ll2004l ) 
suggesting the presence of obscuring matter along our line of 
sight to the BLR, possibly the putative torus. We therefore 
observed this galaxy as a comparison for the Seyfert Is, but 
without the complication of broad emission components on 
both t he permitted l ines and the FHILs. 

As lErkens et al.l lfl997i ) had already found, this Seyfert 
2 contains strong FHILs. We chose to observe this galaxy 
at higher spectral resolution and S/N levels, to determine 
if these lines had any structure in their profiles that might 
indicate outflows or turbulence in the FHIL emitting regions 
of Seyfert 2s. 

The various line profiles extracted from the spectrum 
of Mrk573, with their associated Gaussian fits are shown in 
the last panel of the all the figures discussed previously. 

3.4-1 Narrow Emission Lines 

As expected for a Seyfert 2 spectrum all the emission lines 
in Mrk573 are considerably narrower than the permitted 
lines observed in Seyfert 1 galaxies. However the strong 
[O III] and Balmer emission lines show signs of symmet- 
ric, broad wings and were modelled with multiple Gaussian 
components, while all the weaker lines were satisfactorily 
fit with a single component. Excluding the broad wings of 
the strongest lines, all the emission lines had similar veloc- 
ity widths of FWHM ~300km/s. We also detect strong He 
IIA4686 emission in this galaxy. 

3.4- S Forbidden High lonisation Lines 

Of note is the large number of FHILs that we detect in 
Mkn573. They are comparable to the most highly ionised 
Seyfert 1 spectrum in our study, Ark564. In addition to [Fe 
VII]A6087 both [Fe VII]A5159 and [Fe VII]A5722 are well 
detected. In contrast to the trend found in the Seyfert Is, 
the highest ionisation lines in Mkn573 show no sign of asym- 
metry or wing broadening, and have thus been fitted with 
a single Gaussian. The resulting profiles of [Fe X] and [Fe 
XI] have widths and velocity shifts similar to that of [Fe 
VII]. We note that these lines appear to be redshifted rel- 
ative to [S II], but are consistent with the velocity of the 
Mgb absorption lines. We do not detect [Fe XIV]A5303 in 
this object. 



4 DISCUSSION 

As described in detail in the previous sections, we have been 
able to decompose many of the optical emission lines into 
multiple components that span a range in velocity widths 
and velocity shifts- indicating that the spectra of Seyfert 
galaxies cannot be considered simply in terms of narrow 
and broad emission lines. This is clearly shown in Fig. (2] 
in which we display the range of the component widths for 
each line species which we will discuss in detail. 

While traditionally the permitted lines of Seyfert galax- 
ies were considered to be a blend of emission from both the 



broad and narrow line regions, we find evidence of a very 
broad (2500km/s < FWHM < 7000km/s) third component 
in 6/9 of the Seyfert Is (all six of them NLSls) we observed. 
This extra component w as previously noted in the NLSl 
galaxy KUG1031-I-391 bv iMason et all l|l996l) with a mea- 
sured FWHM (~2500km/s) more t ypical of the broad lines 
in a Seyfert 1. It was suggested bv lGoncalves et al.l (|l999l ) 
that the broad wings of the Balmer lines in KUG1031-(-391 
were better fi t with a single Lorentzian profile, although 
iDietrich et all (|2005l ) have since shown that this function's 
profile overestimates the core of the line, while underestimat- 
ing the flux in the wings, in their study of 12 NLSl galaxies. 
In our appendix, we demonstrate a 3 Gaussian component fit 
is superior to a Lorentzian-plus-Gaussian fit. These models 
reveal a very broad component in a further 6 NLSl galax- 
ies, in which we measure broad components of similar widths 
when fitting both the Ha and H/3 lines independently. We 
note that care must be taken when interpreting lines with 
multiple components. It is unlikely that the dense emitting 
gas around the central engine can be simply divided into 
discrete velocity- width bins (in this case, three). The true 
situation is that we see an amalgam of emission from a en- 
semble of clouds along our particular line of s ight with a 
wide range of velocity widths (see, for example, lOsterbrockl 
Il993h . 

As discussed bv IDietrich et all (|2005h . the fact that we 
see emission line components with widths similar to those 
of the broad emission lines of BLSl galaxies indicates that 
emission from this BLR may be reduced in NLSls, rather 
than absent altogether. This could imply either: 

• that we should no longer assume that black hole masses 
in NLSls are smaller than in BLSls. This may appear to 
contradict other evidence of reduced black hole masses in 
NLSls, including values derived from stellar velocity disper- 
sion (cr) and X-ray variability measurements. It has been 
claimed, however, th at the validi t y of t he Mbh-o relation 
shown in Figure 2 of iBarth et al.l f2005) is uncertain at low 
values of values of a. This is particularly evident when one 
considers that the Mbh estimate for M33 is a hard upper 
limit (and more li kely that this ga l axy d oes not contain a 
central black hole; iGebhardt et al.l l200lf ) and the mass of 
the black hole in NG C4395 has been revis ed upwards by an 
order of magnitude (jPeterson et al.l [20051 ). Taking this into 
account, it is possible that the Mbh-ct relationship shows 
an upturn at low values of a, strengthening the claim that 
we may expect broader components in the Balmer lines of 
NLSls. Also, although detailed long term X-ray variabil- 
ity measurements can be used to derive bl ack hole masses 
using their power density spectra (see e.g. iMcHardv et al.l 
(2000)), isolated short term X-ray variability may not be a 
reliable indicat or of the black hole mass. This has recently 
been shown bv lReeves et al.l l|2002l ). 

• that NLSl galaxies do contain lower mass black holes, 
but the radius of the broad line region scales in proportion 
to the mass of the black hole. This would then mean that 
other factors determine the observed widths of the lines. 

For both these explanations, we would need a theory 
to explain the reduced relative luminosity of the very broad 
component. This could be disruption by the inner parts of 
an energetic wind originating from within the accretion disc 
that would have the largest effect on a nearby broad line 
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Figure 2. The component widths of each of the hnes discussed in detail in this paper. The vertical height of each block represents the 
extent of both the systematic and random errors associated with each line model component. The shaded areas represent the regions 
between the 'narrow' and 'intermediate' and 'intermediate' and 'broad' widths, averaged between the Ha and H/3 components (only one 
shaded region is shown in cases where the permitted line fit required only 2 components). This clearly shows that the FHILs tend to 
have widths between those of the traditional broad and narrow emission lines. 



region, or the strong UV/soft X-ray flux seen in the SED of 
many NLSls in the form of the soft-excess, or a mixture of 
both. 

We find that the permitted Hehum line profiles differ 
considerably from the Balmer lines for the sub-set of galax- 
ies (6/7) in which the H elium line can be measured. This 
fact has been reported bv lShudeil ()l982l ). and has since been 
briefly ment ioned in the course a number of permitted line 
studies (e.g. lAlmog fc Netzeij[l989l . iDietrich et al.ll2005h . in 
which the broad component of the He I A5876 line was found 
to be broader than the corresponding Balmer line compo- 
nent. We too find that this is the case in three of the galax- 
ies in which this line was observed. This suggests some form 
of segregation of the broad line region, in which different 
species of permitted lines are emitted preferentially in re- 
gions with differing conditions and kinematics. This is also 
suggested by the result s of the locally , optim ally emitting 
cloud (LOG) models of [Baldwin et al.l l|l995l ). Because we 
are uncertain of the processes that cause broadening of the 
permitted lines we cannot rule out diverse explanations, for 



example, that a greater proportion of the He I emitting re- 
gion lies closer to the central engine than the Balmer line 
emitting region, because of the higher ionisation of Helium 
compared with Hydrogen, or alternatively whether the He- 
lium clouds are undergoing greater turbulent motion or are 
outflowing/infiowin g at high e r spe eds as has been discussed 
by e.g. CoUin-Sou ffrin et al.l (|l988l ). However, caution must 
be taken when interpreting the widths of the He lines as 
we have shown that the level of the underlying continuum 
can have a considerable effect on the properties of the fitted 
components. 

Due to the difficulties with blending, we have been un- 
able to measure the broad component of the He HA4686 line 
to the same degree of accuracy as He I and Balmer lines. 
Despite this, there is some evidence that it has a broad 
component with a FWHM which is significantly less than 
the corresponding He I and Balmer line components (see, 
in particular, Ark564, Mrk335, NGC7469 and VHZW118 in 
Fig. [2]). If this is the case, and is not the result of measure- 
ment errors or Fe II removal, then this runs counter to the 
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Figure 3. Comparing tlie profiles and velocity shifts of tlie high ionisation lines. We include the narrow [O III] line as a comparison. 
Profiles shown are (from bottom) [O III]A5007, [Fe VII]A6087, [Fe X]A6374 and [Fe XI]A7892. The dotted fines indicate removed features, 
such as the [O IJA6363 line that is often blended with the [Fe XJA6374 line and the blue wing of the Ha line. It is generally seen that 
line widths and shifts increase with ionisation potential. 



predic tion from photoionisation models of iKorista fc GoadI 
l|2004l - that the broad components of the permitted lines 
sho uld increase in widt h with increasing ionisation poten- 
tial. iLandt et all l|2007D find a similar situation when com- 
paring the widths of the permitted lines in both optical and 
near-infrared spectra. They speculate that this effect may 
be caused by an accelerating outflow, in which the He"*" gas 
is closer to the source of the ionising continuum than He", 
but is outflowing more slowly as it has yet to be accelerated 
to its terminal velocity. This would be a very important re- 
sult, but clearly much more work is necessary to determine 
the true kinematics of the He I, He H and the Hydrogen 
emitting clouds. 

Evidence of emission components with widths between 
the conventional broad and narrow lines is also seen in the 
forbidden lines of all the objects we observe, in the form 
of a broad blue wing. This is a well known feature that 
has been st udied in depth since i ts first detection and de- 
scription bv lHeckman et aP (|l98lh . in which they attribute 
the blue wings of the [O HI] lines to outflowing gas, which 
is likely to have a density similar to the traditional NLR 



(i.e. rie '-^lO'^cm"^), interspersed with dust that preferen- 
tially obscures the gas outflowing in a direction away from 
us. More recent studies (.Whittle. 1992I. Christopoulou et al.l 



1 19971 . IVeilleux et allbOOll . IVeroii-Cettv et al.ll200ll ) confirm 
that the [O HI] emission is likely to arise from an am- 
bient region at the same redshift and rotational veloc- 
ity as the host galaxy, and an 'outflowing' region that is 
distributed in two wide angle ionisation cones. We note 
that the critical densities of [O HI]A5007, [Fe VII]A6087 
and [Fe X]A6374 (log(ne)^5.8 -7.6, ~9.7, respectively; 
lAppenzeller fc Oestreichenl 19881 ) suggests that the outflow- 
ing gas responsible for [O HI] blue wings is likely to have a 
lower density than the FHIL outflows. 

Strong evidence of gas with kinematics between the 
classic narrow and broad emitting re gions is found in the 
form of the FHILs. As was noted in Erkens et aP (|l997l ) 
and iRodrfguez-Ardfla et al.l l|2006l ). we find that as well as 
being significantly broader than the lower ionisation forbid- 
den lines, these lines show increasing blueshift with increas- 
ing ionisation potential. This is clearly shown in Fig. [S] in 
which we plot the line profiles in terms of velocity shift from 
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systemic for the [O III]A5007 and the Fe FHILs, which sug- 
gests an outflow that is possibly related to that responsible 
for the broad, blue wing of the [O III] lines. 

We note that, where it is measured, the width of the 
[FeVII]A6087 line is generally between that of the narrow 
and intermediate Balmer line components. This, combined 
with findings that this line also tends to be blueshifted with 
respect to the permitted and low ionisation forbidden lines 
suggest this line is emitted from a kinematically distinct re- 
gion. We find that the relative intensity of the [Fe VII]A6087 
is reasonably constant among our small sample, at 2-4% fiux 
of the H/3 line in all the galaxies we observed. A larger, more 
homogenous sample would be required to fully determine 
correlations between the [Fe VII] line properties and AGN 
classes, such as NLSls and BLSls. 

Our fits to the [Fe X]A6374 line suggest that this species 
is formed in yet another kinematic a lly dis tinct region, as was 
noted by iRodn'guez-Ardila et al.l l|2006l ) in nearby Seyfert 
1 and Seyfert 2 galaxies. Although approximately twice as 
strong as [Fe VII]A6087 in all the galaxies in which it is de- 
tected, this line is more poorly defined because of its much 
broader profile and its blending with other emission lines. 
In many of the galaxies it has a broad component that ap- 
proaches the FWHM of typical of broad permitted lines, and 
shows little or no evidence of a narrow component, possibly 
implying that the Fe®^ ion is only created in regions near to 
the central source of ionisation, which woul d agree with the 
predic tions of the photoionisation models of iFerguson et all 
|l993)- As found for [Fe VII] the relative intensity of the [Fe 
X] line appears to be similar in all the galaxies we observed 
(~8% H/3, or 2x[Fe VII]A6087). We note that the almost 
constant value of [Fe X] A6374/[Fe VI I ]A608 7~2 appears at 
odds with the findings of lNagao et al.l l|2000l) . who generally 
found lower values for this ratio in a larger sample of AGN 
(including NLSls, BLSls and Sy2s). We note, however, that 
the systematic errors associated with the fitted continuum 
level has a significant effect on the measured intensity of the 
broad components of these weak lines, which could account 
for these differences. 

In contrast to generally similar [Fe X]A6374 relative line 
intensities among our sample, we find that the [Fe XI] /[Fe 
X] ratio covers a wide range e.g. from ~1 in 1111934-063 to < 
0.03 in NGC7469 (taking 90% upper confidence limits for the 
[Fe XI] line). This is despite the fact that these two species 
have similar ionisation potentials (235eV and 262eV for [Fe 
X] and [Fe XI], respectively). The origin of this wide range 
of ratios is not well understood. How ever, it is in agr e ement 
with the photoionisation models of IFerguson et al.l ()l997l ) 
who point out that the [Fe XI]A7892 line is highly sensitive 
to the density and ionisation parameter of the FHIL emitting 
region. In addition, if we assume t hat [Fe XI] is produc ed 
by photoionisation (as suggested bv lPenston et all 1 19841 ') it 
is possible that the range in the line ratio is also a function 
of the continuum shape in this energy region. 

Despite the difficulties in extracting the [Fe XIV]A5303 
mentioned in i|3.3l we are able to qualitatively confirm the 
presence of this line in two of our sample, Ark564 and 
1H1934-063. These galaxies also display strong [Fe X]A6374 
and [Fe XI]A7892 lines, and we are able to show that the 
blueshift of the [Fe XIV]A5303 line is similar to the blueshift 
of these other lines, giving us further confidence in this de- 
tection. However, it should be emphasised that the system- 



atic errors associated with this measurement caused by the 
removal of the Fell blend is unknown and are likely to be 
large, we therefore make no attempt to interpret the inten- 
sity or width of this line. 

Finally, turning to the case of the only Seyfert 2 in our 
sample, it is clear from Fig.[2]that the highest ionisation lines 
in Mrk573 trace a kinematically different region than both 
the permitted and other forbidden lines. Apart from a weak, 
broad component of the [Fe VII]A6087 line, the widths of 
FHILs are all relatively l ow (FWHM 400km/s). The se re- 
sults agree with those o f lAlmudena Prieto et al.l (|2005| ) and 
iRodrfguez-Ardila et al.l (|2006l ). who found that there is sig- 
nificant [Si VII] (a FHIL observed in the near-IR) emission in 
regions up to ~150pc from the nuclei in four nearby Seyfert 
2 galaxies. However, it is difficult to reconcile their results 
with our findings for BLSls and NLSls in terms of the stan- 
dard AGN model. Strong [Fe XI]A 7892 is detected in all 
four S eyfert 2 galaxies observed by iRodrfguez-Ardila et al.l 
(j2006h and in our work. However, this line is detected in 
less than half of the Seyfert Is that have been observed by 
us and in lErkens et all (|l997l ). If, using the standard model 
we expect that the NLRs of both Seyfert 1 and Seyfert 2 
galaxies are essentially the same then we would expect from 
the Seyfert 2 results to detect a narrow [Fe XI]A7892 line in 
a much larger proportion of Seyfert 1 galaxies. Of course, 
sample selection effects may be important and it will be 
necessary to observe larger samples of the various classes 
of Seyferts at comparable signal to noise to make further 
progress in this area. 



5 SUMMARY 

The main conclusions from our study can be summarised as 
follows; 

• The emission lines in the spectra of both BLSls and 
NLSls appear to be emitted from more zones than the tra- 
ditional narrow and broad line regions. This is particularly 
evident in the profiles of the FHILs, but also in the profiles of 
the broad, permitted lines of Ha and H/3. The lines from the 
FHIL region typically have a FWHM of 1000 - 3000km/s, 
and are significantly narrower than the broadest components 
of the permitted lines observed in the same object. 

• All of the NLSls in our sample show evidence of a very 
broad component in their permitted lines. This component 
has a FWHM that is more typical of that observed in the 
broad lines of BLSls, i.e. ~3000km/s. Affhough the FWHM 
measurements are a product of our line modelling process, 
the point that the broad wings of the permitted lines ex- 
tend to up to 3000km/s in even the 'narrowest' NLSls in- 
dicates that there exists emitting clouds with much higher 
relative velocities than has typically been assumed in this 
sub-set of Seyferts. However, this component appears to be 
greatly suppressed in the case of NLSls. We speculate this 
suppression might be either a result of the destruction of 
the standard BLR by an energetic wind or by the increased 
high energy photon fiux of NLSls, or both. We propose two 
methods of production of this extra broad component: 1) a 
higher mass black hole than is expected for NLSls and a 
BLR radius more typical of BLSls, or 2) a low mass black 
hole and an equivalently smaller BLR radius. 
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• In a number of our sample the He I A5876 line con- 
tains a component that is broader than the broadest Ha 
and H/3 components. This may be an indication that the 
permitted He lines are preferentially emitted closer to the 
central engine than the H permitted lines, or that a differ- 
ent mechanism is causing the broadening of these emission 
lines. 

• We detect [Fe XIJA7892 in a few cases where the data 
has sufficiently high signal to noise. Despite [Fe X] and [Fe 
XI] having similar ionisation potentials, the [Fe X]/[Fe XI] 
ratio covers a wide range within our sample. It is therefore 
possible that this ratio is very sensitive to the physical con- 
ditions of the emitting region as well as the shape of the 
ionising continuum. This relationship will be investigated in 
a follow up study. 

• Strong, narrow FHILs are^ also detected in the Seyfert 
2 Mrk573 (as also noted by Rodriguez- Ardila ot al. 200^. 
Since it has been shown that this AGN contains a hidden 
BLR, this result indicates that FHILs can be formed in re- 
gions beyond the extent of the putative dusty molecular 
torus. 
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APPENDIX A: MULTIPLE COMPONENT FITS 
TO THE PERMITTED EMISSION LINES 

A statistical goodness-of-fit test is required in order to 
conclusively determine whether multiple Gaussians or a 
Lorentzian plus Gaussian model best represent the shape 
of the observed broad permitted emission line profiles. 

In the literature the principal statistical test for this is 
the fit statistic, which takes account of both the data 
points and their associated errors. When we apply this 
statistic to our Balmer line data we find that the values 
of the three-Gaussian (3G) models are significantly lower 
than those of the Lorentzian-plus-1- Gaussian (L-I-IG) mod- 
els. For example, in the case of the H/3 line in Ark564, the x^ 
statistic for the 2 Gaussian is 6077.97, L-I-IG: 2249.38 and 
3G: 2186.66. The best-fitting L-I-IG and 3G models both 
provide a much better description of the data than the two- 
Gaussian model. Comparing these two models, we find that 
the 3G model improves the fit by Ax^ = 62.72 with the 
addition of just three additional free parameters. The prob- 
ability of three extra parameters improving the x^ by as 
much as 21.10 purely by chance is 0.01%. The Ax^ is at 
least this large for each line in tables 2-4 that is described 
with three components. Thus, we assert that the 3G model 
provides a better fit with >99.99% confidence. 

We emphasise that the reduced x^ of all of the above 
models is S> 1, so none of these provide a formally good fit 
to the data. This is an inherent problem when fitting simple 
models to large-scale features in high signal-to-noise data. 
The random error associated with the flux in each spec- 
tral bin is small and thus the potential contribution of any 
individual bin to the total x^ is large. Small 'bumps-and- 
wiggles' in the data that are not reproduced in the model 



will thus contribute strongly to the total x^ ■ Although each 
of these deviations may well be a real feature of the source, 
we are not presently interested in complicated models that 
could reproduce such details. The inevitable consequence of 
leaving such features unmodelled is that the reduced x^ of 
the best-fitting model will invariably be larger than 1. In 
this paper we are only interested in the large-scale struc- 
ture of the lines. If we want our relatively simple models to 
fit the spectra well we need to minimise the impact of the 
small-scale bumps-and-wiggles. This can be accomplished by 
smoothing the data. So long as the smoothing radius is not 
too large, the gross features in which we are interested will 
be preserved. 

To further test the L-f IG and the 3G models we have 
applied both to a sequence of increasingly smoothed spectra. 
We convolve each spectrum with a range of Gaussian kernels, 
varying the Gaussian FWHM from l.oA to 25. OA in 0.5A 
steps. To the lines in each we fit the L-I-IG and 3G models 
and measure the resulting reduced x^ ■ We note that by using 
the reduced x^ ^6 take into account the greater number of 
free parameters (and therefore greater fiexibility) associated 
with the 3G model. As the smoothing radius increases, (1) 
the small-scale structures in the data get washed out, (2) 
the information content of the spectrum decreases, and (3) 
the reduced x^ of the best-fitting large-scale line models 
decreases. 

With the appropriate amount of smoothing, the reduced 
X^ of a model will approach 1.0 and the model can then be 
said to provide a good description of the (remaining) data. 
Any additional smoothing will decrease the reduced x^ be- 
low 1. When this happens the model is no longer appropriate 
because the spectrum has lost so much information that it 
no longer requires a model of this complexity. In effect, the 
model has too many degrees of freedom for such a simplified 
spectrum and its use over-interprets what is left of the data. 
We find that for the lines measured, the 3G models require 
less smoothing than the L-I-IG models to reach a reduced x^ 
of 1. We therefore consider the 3G model to be the more ap- 
propriate model because it describes the spectra well when 
less of the original information content has been lost. 

As an example, fig. lAll shows a comparison of the 
reduced-x^ values for both the 3G and L-I-IG models ap- 
plied to the H/3 profile of Ark564 as a function of smooth- 
ing width. The reduced-x'^ of the 3G model reaches 1 with 
a smoothing kernel of ~4.0A FWHM whereas that of the 
L-I-IG model requires a kernel of ~5.6A FWHM. Thus the 
3G model refiects more of the information contained in the 
original spectrum. Quantitatively, when the 3G model fit- 
ted to the smoothed spectrum has a reduced-x^ of 1.00 the 
reduced-x^ of the L-I-IG model is 1.35, indicating that the 
former is a better representation of the data with >99% con- 
fidence. Moreover, the fact that the reduced-x^ curve of the 
3G model is always below that of the L-I-IG model indicates 
that the 3G model always provides a better representation 
of the data regardless of the level of smoothing applied. 

It is important to note, however, that care must be 
taken to choose an appropriate level of smoothing; i.e. 
smoothing such that the weak, narrow, 'bumps-and-wiggles' 
are removed from the data without strongly affecting the 
broad line profile itself. Our primary concern is to ensure 
that the smoothing does not significantly alter the narrow 
component that makes up the core of the broad Balmer line 
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Figure Al. The effect of smoothing on the reduced-x^ statis- 
tic for the Ark564 H/3 Hne (here, smoothing length refers to the 
FWHM of the Gaussian smoothing kernel). The dotted line cor- 
responds to the L-I-IG fit, and the solid line refers to the 3G fit. 
It is evident that, irrespective of the smoothing length used, the 
3G fit is statistically better than the L-I-IG fit. That the reduced- 
X^=1.0 point is reached by both fitting methods at smoothing 
lengths significantly below the maximum smoothing length re- 
vealed in I A2I confirms this is still within the regime at which the 
narrow component of the H/3 remains in tact after smoothing 
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Figure A2. Change in x'^ after the removal of the narrow com- 
ponent of the H/3 line of Ark564 when fitted with a 3G model. 
The dotted line shows the 90% confidence level that the narrow 
component is detectable when performing the fit. Using smooth- 
ing lengths significantly less than point at which the two lines 
intersect effectively ensures that we are smoothing out the small 
deviations discussed in the text, but retaining the overall profile 
of the H/3 line. 



upper limit on the appropriate smoothing value for an indi- 
vidual line profile depending on how stringently we wish the 
detection of the narrowest component to be. Fig. IA2I shows 
the change in when the narrowest component is removed, 
as a function of smoothing width, in the case of Ark564. For 
each emission line we limit our chosen smoothing scale to 
be much less, that is by at least a factor of 2, than the 
maximum smoothing scale at which the narrow component 
can still be detected at the 90% confidence level. This is to 
ensure that the narrowest component of the line is not sig- 
nificantly affected by the smoothing process. This places an 
upper limit on the appropriate smoothing width of ~4.75A. 
We note that the 'best-fit' smoothing width for Ark564 is 
^4A, which is comfortably below this upper limit. 



profile. To ensure that we were still able to measure the prop- 
erties of the narrowest component after smoothing, we fit a 
3G model to the unsmoothed data, recording the properties 
of the narrowest component. Again, we then convolve the 
line profile with a Gaussian kernel whose width increases 
from I.OA FWHM to 25.0A FWHM, in steps of O.SA. At 
each step the smoothed data are re-fit with a 3G model in 
which the width of the narrowest component was fixed at 
the width of the narrowest component, based on the fit to 
the unsmoothed data. After calculating the we then sub- 
tracted the narrowest Gaussian from the model. The change 
in that results from the removal of the narrowest compo- 
nent of the model reveals how well the narrow component 
is detected for a given level of smoothing, and so places an 
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Figure A3. Individual line plots and component fits. For this and all the following plots of the line profiles, we only show the spectra 
for the galaxies in which the line is detected and measured. In each panel, the vertical dashed lines show the rest wavelength of the line 
of interest. Dotted lines reveal where we have removed a blended line to assist component fitting. Colour coding- Red: Full line model. 
Magenta, Green, Blue: Individual lines components (broadest to narrowest). Orange: Blended component not associated with measured 
line (e.g. [N II] blended with Ha). For the case of H«, shown here, we have removed the [N II] lines in Mrk573, as described in the text. 
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Figure A4. H/3 & [O III]- Dotted lines show where we have removed the [O IIIJA4959 Hne in order to more accurately fit the H/3 lines. 
In the case of VIIZW118, we have also removed the He IA4922 that forms a strong red shelf in this spectrum. 
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Figure A6. He IIA4686- We have deblended this line from the unknown line at ~4670A, which is clearly a separate feature in a number 
of our targets (e.g. IIZW136, VIIZW118) 
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Figure A7. 
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Figure A8. [Fe XJA6374- Dotted lines here are used to show where we have removed both the [O I] hne that is blended with the [Fe X] 
lino. Rodward of some of the [Fe X] lines we have removed the blue wing of the Ha line that contaminated the [Fe X] line in some of the 
'broader' NLSls in our sample. 
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Figure A9. [Fc XI]A7892 
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Figure AlO. [Fe XIV]A5303- For Ark564, we have attempted to deblend this hne from the [Ca V] hne at 5309A (shown here in orange) 




